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V-3. 


Introduzione 


Due anni fa,nell'ambito di questi seminari, Graffi [4] ha mostrato co 
me l'equazione di Schrédinger stazionaria può essere ricondotta ad un'equazio- 
ne di Hamilton-Jacobi più correzioni in f ed ha dimostrato come, nel caso del- 
l'operatore di Schrédinger dell'oscillatore armonico non risonante perturbato 
con un potenziale polinomiale, il k-esimo coefficiente dello sviluppo di Ray- 
leigh-Schròdinger di ciascun autovalore si esprima come somma del k-esimo ter- 
mine della serie normale di Birkhoff più correzioni in fi; questa rappresenta- 
zione degli autovalori è ovviamente molto interessante per trattare il proble 
ma del limite semiclassico (v. [2]). 

Vogliamo qui riprendere il discorso mostrando una nuova applicazio- 
ne di quel metodo al caso dei rotatore forzato con una forza periodica, proble 
ma già trattato in [3]. Rispetto al caso esaminato da Graffi [4], qui si pre- 
senta una difficoltà ulteriore nel caso in cui la frequenza propria del siste- 
ma è risonante con quella del termine forzante. Alla trattazione del problema 
premettiamo quindi un breve resoconto della teoria canonica delle perturbazio 


ni nel caso risonante. 


Teoria delle perturbazioni classica 


E' ben noto che una hamiltoniana H:W + R, definita in un sottoin- 
sieme W del fibrato cotangente T*5 ad una 2-varietà 2, si dice (canonicamente) 
integrabile se esiste una trasformazione canonica C:W + vxt* tale che 


1 


(1) H(C (A,9)) = h(A) 


x i É <a -, ; : 
ove V è un aperto di R' e T è il toro unidimensionale, mentre 1a trasformazione 


A 
C si dice canonica se conserva la 2-forma simplettica w = :D UP bd in TS, 
J5l 


In questo caso le equazioni di Hamilton per la nuova hamiltoniana 


definita da (1) sono 


è = 2A) 
À=0 i 
#(0) = è, 
A(0) = A; 


esse si integrano immediatamente 


d(t) w(AT + do 


A(t) = A 


(o) 


“ 


avendo posto w(A) = di (A;)* 11 moto (2) risulta così essere quasi periodico. 
La teoria classica delle perturbazioni tratta il problema del moto 


generato da una perturbazione di una hamiltoniana integrabile 


(3) H_(a,p) 3 H(q,p) + eV(q,p) 


In generale la trasformazione canonica C non integra la hamiltonia- 


na (3), infatti 
(4) H_(C"°(A,9)) = h(A) + ef(A,6) = h(A,6). 


E' stato inoltre provato da Poincaré [8; vol. 1 ch. V] che, in gene- 
rale, è impossibile trovare una nuova trasformazione canonica Co analitica in 
e per e in un intorno dell'origine, che trasformi la hamiltoniana (4) in una di 


pendente solo dalle azioni A. 


(5) h_(C-2(AL,8'_)) = (81). 


E € 


In generale però per ogni kE N è possibile trovare una trasformazio 


ne canonica Cc. tale che (5) valga fino all'ordine eb nel senso che 
-l ' ' HI ' k 
(6) h_(0- (A',9')) = h_(A')to(e°) per ec +0. 


Questo può essere ottenuto risolvendo le equazioni di Hamilton-Jacobi per la 
funzione generatrice M_ all'ordine ALTA II ch. IX]cfr. anche [1. 8 3]). 


94 


(7) ho 39) = RA") +ole 


M 


Per questo si scrive la funzione generatrice come un polinomio di 


grado k in e 

(8) H_(A',0) = A'g + ell, (A'.9)+...+ e" Wi, (A',9); 
si cerca poi la nuova hamiltoniana hi nella forma 

(9) R_(A") = N_(A") tell (A') + ...4 e"W.(A") 


ed infine si sviluppano entrambi i membri di (7) secondo le potenze di e otte- 


nendo 
(10) hCA') + e Pala!) Pa ef... def(A',g)tel... = 
= N CA") + eN (A') + "a 


Uguagliando ora i coefficienti delle successive potenze di e în (10) si otten- 


gono le equazioni 
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E 
-_ 
ba 
_ 
Qu 
e 
" 


N_(A') - F(A',6) 


w(A') —I N,(A") + P;(A',0) $ © Bk 


dW 
essendo P.(A',0) un polinomio in LO, con h = 1,...,j-1. 


I] caso non risonante 


Le (11) si risolvono per ricorrenza su j utilizzando gli sviluppi 
di Fourier; in particolare per la seconda si cerca lo sviluppo di Fourier del- 


la funzione incognita 


(12) Wi (A'6) = DO my p(8) gi”9 
ne" 


in funzione dello sviluppo di Fourier del termine noto 


(13) £(A'19)= )_ f_(A') E 


A n 
nez 


Inserendo (12) e (13) in (11) si ottiene 


(14) i<w(A').n> w, (A') = -f(A') n7#0 


La n 


1 
(15) N_(A') = f_(A') = f(A',9)dò 
ì d (27)? J, 


V-7. 
Evidentemente (14) può essere risolta solo se 
(16) <w(A'),n> #0 Yn #0; 


in tal caso essa dirà luogo ad una serie convergente in (13) se w(A') soddi- 


sfa una condizione diofantina del tipo 
(17) <w(A'),n> > C [n|T® n#0 


essendo C>0 ed a>0 due costanti indipendenti da n. 
Chiaramente (16) è la condizione che le frequenze w(A') non siano 
risonanti mentre (17) può essere interpretata come una condizione di lontanan- 


za dalle risonanze. 


Il caso risonante 

Nel caso risonante,in cui la condizione (16) non è soddisfatta, è 
evidentemente impossibile risolvere (14) ed è quindi impossibile trovare una 
funzione generatrice della forma (8) tale che la trasformazione canonica da 
essa generata fornisca una hamiltoniana della forma (6)-(9). Anche in questo 
caso tuttavia è possibile trovare una trasformazione canonica per la quale val 
ga ancora (6) e questo si ottiene applicando i1 seguente metodo dovuto a Boh- 
lin e Poincaré [8 vol. II ch. XXI]. 


Supponiamo dunque che esista mez” sì m#a0 


<u(A'),n>#0 Wn#jm Wjez, 


V-8. 
e cerchiamo una funzione generatrice della forma 
' ' 3 ' L] 3/2 4 
(19) M_(A',9) = <A",99 + e°W(A :9)+ eWp(A',9) + e W3(A',0) + 
k 
+ E Max (A 39). 


Allora è possibile determinare le funzioni Mircea e NoN: 3 N, in modo 
che la trasformazione canonica C_ generata dalla funzione M_ (19) trasformi 
la hamiltoniana h_ (4) nella forma (6),dove la nuova hamiltoniana h data anco 
ra dallo sviluppo (9). 

Inserendo (19) e (9) in (7) e sviluppando secondo le potenze di Je 


si ottiene 


aW aW dW, WU 
1 Ù ' 
papa i ij da ay 10 HA) _L_L + 


9A 90 9A 99 2 aA3A 96 94 
(20) 
+ 8/2 2b(q1) 23, 3/2 2° (n°) pel ARE 
9A dò 9A9A dd 96 
3/2 9df 2 


OM, 2 
0. + ef(A',9) + e AA ,9) o €... 3 NSA )te N (A te... 


Questa è l'analoga di (10) e permette di ottenere le analoghe di (11) 


h(A') = NCAA ) 
dW 
w(A') — = 0 
99 
9U 2 9W, 9W 
(21) ata DDR) LA Lula . 
w(A') 20 +2 3A 3A dé dd f(A',9) = N,(A ) 
Mii FCI ho Lap SENT PT 
di dd 2 9A 9A 29 39 3A 385 sé DA *9% 34 
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La prima equazione delle. (21) coincide con la prima delle (11), men 
tre a differenza della seconda equazione delle (11) la seconda delle (21) non 


determina la funzione Mi: da essa è possibile solamente dedurre che 


(22) W(A',6) = H_(A5<m,6>) 


ove la funzione N (A',0) è periodica di periodo 2r in 1. Per determinare tale 


funzione si ricorre quindi alla terza delle (21) ponendo 


(23) W_(A',6) = W_(A',4) + N_(A',<m,4>) 

ove 

(24) N_(A',9) = Ù, My nAVe!" 
nez* 


nè jmVjezZ\{0} 


(25) U_(A',<m,0) = LG, (A') elim 
jezo} 


Sostituendo (25) e (24) in (23) e quest'ultima nella terza equazio- 
ne della (21) si vede subito che la funzione LA rimane indeterminata mentre è 
possibile determinare N,» N, ed N, separando le frequenze risonanti da quelle 


non risonanti; per questo poniamo 


(26) î(A',9) = a f_(A') pe 
nez* 
n#jm, j € Z\{0} 


(27) tas) = Lor, (A') eli 


jentor I" 
onde 
(28) lA") 24 (At) = RCA) 
e 
(29) 12° (A') ne 23, #(A',0) = Î(A") 


E' evidente che l'equazione (28) determina le funzioni Wp ed N, 


nello stesso modo in cui la seconda equazione di (11) determina HW, ed Nj5 le 


formule che forniscono N, ed Î. sono del tutto analoghe a (14) e (15). Per 


1 
quanto riguarda invece l'equazione (29), posto t = <m,é> e 


B(a)-L AIA) min > 
otteniamo 
dii, 5 
(30) B(A') (FI) + HA',t) = N (A'). 


Si cerca ora di determinare N, (A) in modo che l'equazione non lineare (30) ab- 
bia una soluzione N» periodica di periodo 2r nella variabile 1 per ogni A'. 
A questo punto è possibile risolvere le rimanenti equazioni ricorsi 


vamente: la quarta ad esempio fornirà Ù, e DA etc. 
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Equazione di Schròdinger ed 
equazione di Hamilton-Jacobi 


Consideriamo ora l'equazione di Schridinger 
2 1,2 8° 2,2 
(31) if rt dalla Dai p3vte Va,wt)y = H_(t)v wveL'(T°) 


da 


dove il potenziale ha la forma 


(32) V(a,8) = V(a) COSB » 
con 

eo : 
(33) V(a) =). vie"; 


supporremo che la funzione (33) abbia un prolungamento analitico in un intorno 
complesso del toro T nel senso che esiste una costante E>O0 tale che la serie di 


Laurent 


(34) È. È Vv” 


converga nell'intorno 


Ce {zec: et < |z]< a5} 
Evidentemente H_(t) è la hamiltoniana del rotatore perturbato dal termine for- 
zante periodico (32). 


E' noto il seguente risultato di VYajima. 


Lemma [9]. La funzione w è una soluzione dell'equazione d' evoluzio- 


ne (31) se e solo se la funzione 
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è soluzione dell'equazione stazionaria 


2,2 


hi 8° sla a 
a + eV(a,B) - ifiw PAL) E$ » 


(35) A 


cioè $ è un'autofunzione dell'operatore 


2 


(36) raLbd 


2 
- d 
i — SOI — + 
5 2 2 iNw 98 eV(a,8) 


da 


relativa all'autovalore E. 
11 legame fra l'equazione (35) e la teoria delle perturbazioni clas 
sica è fornito dal seguente risultato provato da Graffi e Paul [2] e da Graffi, 


Paul e Silverstone [3]. 
Lemma. Sia 


(37) Siae) att) 


Allora $ è un'autofunzione dell'operatore T, con autovalore E 
(38) T_9 3 Ed 


se e solo se la funzione W è soluzione dell'equazione 


“a: 
1,34,2 36 iti aW_ 
(39) 250) tg * eV(a,B)- 7 È 
da 
ih aÉu 
E' subito chiaro che, a meno del termine - "n (39) è l'equazio 
da 


ne di Hamilton-Jacobi per la hamiltoniana 


(40) T(AB30,8) = 3 AÎ + uB + eV(a,8) (A,B,a,8) € TeTÒ, 
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questa si presenta proprio nella forma (4), cioè come perturbazione della ha- 
miltoniana integrabile SAR già espressa in variabili azione-angolo, con il 
potenziale V(a,8) (32). La hamiltoniana (40) è da un lato la controparte del- 
l'operatore di Schròdinger T. (36) e dall'altro è la rappresentazione indi- 
pendente dal tempo della hamiltoniana del rotatore sottoposto ad una forza pe 
riodica 


(41) X(A,a,t) = 3 A° + eV(a,ut); 


quest'ultima corrisponde evidentemente all'operatore di Schròdinger dipendente 
dal tempo H_(t) (31). 


Soluzione dell'equazione di Schrédinger 
mediante la teoria delle perturbazioni 


Ci proponiamo ora di trovare una soluzione approssimata dell'equa- 
zione (39); per (37) e (38) essa fornisce le autofunzioni e gli autovalori ap- 
prossimati dell'operatore di Schròdinger T, (36). In analogia con (8) e (9) 


cerchiamo una soluzione al primo ordine in e della forma 


(42) E = E, + Ei 
is; 
(43) W(0,8,E,e) i HM (0383E,) + HM (a»B,Ese) 


dove E, e LA sono soluzioni di (39) all'ordine zero in e, mentre E. e W. sono 


1 i 
soluzioni all'ordine uno in e; più precisamente 
2 
1 Wo 2 Mo hi d HW 
(44) raro, ei E 


da 
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e 

aW_ 9W oW 3°u 

o 1 ci ponte i 
(45) a Sa + w se + e V(a,8) - 2 da E, 
da 

La (44) si risolve per separazione di variabili ponendo 
(46) E; b Fal : E,,2 
e 
(47) LACCO E) i Mo,10%E0,1) + Mo, 268»E0,2) 


Questo fornisce subito 


Mo ,16%E0,1) = V2E,.1 a 
(48) 
ad 
Mo ,248»E0,2)  w Egf 


i /W 
“ha a : ; (e e ua . : 
Poiché d'altra parte è necessario che la funzione e sia periodica di perio 


do 2r in a e BR si ottengono le condizioni di quantizzazione 


m 
" 


)? nez 


1 
a (fin 
È = whm mE Z. 
Sostituendo (49) in (48), (47) e (46) si ottiene la soluzione ‘all'ordine zero 
1 2 
(50) E, = (nh) + fim n,mez 


(51) Mo (0383E0) = nha + mh8 
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Per risolvere (45) osserviamo che per (33) il potenziale (32) è svi 
luppabile in serie di Fourier 
+0 


(52) Va,8) = )__ - v_(e' (rete) + gi (ra=8)) 


ov 


Cerchiamo dunque una soluzione di (45) della forma 
- i (ra+8) i (rag 
si (1 (i 
(53) My(@sesEse) =D0, (Wp 1(Esele +, (Gal a 


Inserendo (52) e (53) in (45) si ottengono delle equazioni per i coefficienti 


di Fourier della funzione Hi e la costante E; della forma seguente 


| h_2 De 
(54) i(Ar twt n) ro) Wr,416F>©) "€37 Vi. 
e 
(55). È ce #1 V(a,8)da dg = 0, 


(27) T° 
ove si è posto A = nh. 
Evidentemente (54) fornisce i coefficienti M_ 4] PeP ogni r se e so- 
cdl 


lo se 
Ar+w+ir #0 Vrez; 


in tal caso questi coefficienti risultano di ordine e. 
E' tuttavia subito evidente che il caso di A risonante, in cui esi 
ste meZ tale che Am+ w= 0, è molto diverso da quello non risonante per cui 


Am + w # 0 per ogni meZ. Nel primo caso infatti risulta 


Mail 2% +0 
— fi+0 
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per tutti gli r tali che r = jmcon jEZ. 

Poiché desideriamo ottenere stime uniformi per h>0, è chiaro che 
procedendo neì modo descritto questo è possibile solo nel caso non risonante; 
nel caso risonante è necessario procedere diversamente separando e trattando 
diversamente i modi risonanti e quelli non risonanti nel potenziale (52) e nel 


la funzione incognita (53). 


Il caso non risonante 


Nel caso non risonante le equazioni (54) e (55) insieme alla (53) 


forniscono la soluzione cercata. 


Teorema 1. Siano A,BER tali che esistano due costanti C>0 e y>0 ta 


li che 

-l Y 2 
(56) Ar + us| © < C(Jrl+|s]) (rss) e Z \{(0,0)}. 
Allora per ogni mn eZ ed f > 0 con 


nh 


" 
» 


(57) 
mh = B 


: i ; 2 
esistono un autovalore approssimato E ed una autofunzione approssimata é:T + C, 


nel senso che 
(58) Ta = Ep +R 


ed esiste una costante C>0 tale che se fi è sufficientemente piccolo 


z 
(59) IR(a,6)] sCe ol 
LE (78) Etro) 


Traccia della prova. Abbiamo già mostrato grazie a (46),(47),(48), 
(49),(50)(51) e (55) che una soluzione della forma (37) con Ee W dati da 
(42) e (43) e soddisfacenti (44) e (45) assume la forma 


(60) E=FA° + u 
(61) W(0,8,E,€) = ah + 8BYW(a,8:E,e), 


dove la funzione LA è data da (52), (53) e (54). 
Da (39),(44) e (45) è evidente che il resto R che compare in (58) 


è dato da 
i 2 
(62) R(a,8,E,e) si 2 (7 (038)) d(0,B) 
(0) 
MI 2 
Per pervenire alla stima (59) è dunque sufficiente maggiorare max IT (038) | È 
2 


I 


A tal fine utilizziamo un risultato di Herczynski [5]; precisamente 


h_2,-1 1 ST 
(63) |Ar + w + Dili |a pri lo] 


che permette di ottenere una maggiorazione indipendente da fi; (63) si prova u 
tilizzando la condizione diofantina (56). 
Il caso risonante 


Passiamo ora ad esaminare il caso risonante; fissiamo dunque due co 


stanti A,BER in modo che esista peZ tale che 
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(64) Ap+w= 0 

e scegliamo fjn,m con A= hn ,B=fime 

(65) |Artwt 3 ri |20>0 reZ con r fp. 
Allora vale il seguente risultato . 


Teorema 2. Esistono un autovalore approssimato E ed un'autofunzione 
approssimata é$ : Lai > C della forma 


(66) EE, #E 
e 
(67) $(a,B) = ell? W(a»8) $, (0,8) 
tali che 

LA, = Ept+tR 
con 
(68) 23 Se? c(1+0(1)) max lol I'M, , + 

EST) 12 L'(T) 
2 
+ e Clol , 
ig ng 


dove 0(1) è da intendersi per fi + 0. 
I due addendi di (66) sono 
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dl 22 
(69) E, "2 A+ wB 
e 
(70) E = eV + pf Siti ripe) + Di) per tro 


dove V è il valor assoluto del coefficiente di Fourier ‘o della componente ri- 
sonante del potenziale (52) ed r(p,fi,c) è una funzione limitata; si noti la 
analogia di (69) con (50) e la diversità di (70) da (55). Similmente dei due 
fattori di (67) il primo è del tutto analogo all'autofunzione ottenuta -nel 
caso non risonante mentre il secondo è una funzione di Mathieu la cui espres 
sione è data da (88). 


Traccia della dimostrazione. Abbiamo già osservato come nel caso ri 


sonante (64) non sia possibile risolvere (45) mediante (54) ed ottenere anco- 


ra delle stime che non dipendano da fi; per questo separiamo le componenti non 
risonanti e risonanti pemendo 


(71) V(a,8) = ) LE el(Pot8) % Var e Î(ra+8), 
rez è 
rip 
"p i(pat8) , "p -i(patg) 
(72) Va,8) = # GPIPRTRE 4 va e'*P 
e 
W = i(ra+8) ., -i(rat+g) 
(73) M, (a,8) > L (W_._1 e + La RE e ) 
ràép 


ù ho ij(pat8) 
(74) i (a,g) = L, Wi 05 € 
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Le componenti non risonanti soddisfano un'equazione simile a (45) 


2: DI 

a°A 
o "1 LR 1,4 st 
(75) gd rta È, 


e quindi Ù, ed E, sono ancora dati dalle formule (54) e (55). 


Per quanto riguarda le componenti risonanti richiediamo che 


9H. 9@ au 3°N 3W 
(76) Dei "E +1 (2) 4 Ha,g)= È, 
Poniamo 
(77) pu = "% 
e 
(78) S(u) = N (0,8) 


Osserviamo che la somma dei primi due addendi di (76) è nulla per (64) e se 


V= |v|]ey= argv allora 
| È, Y o%, 

(79) Ya,B)=V cos (2uty) 

e così l'equazione diviene 
1,p as;2 . fi pi 2Îs n 

(80) iP 99 TE + eV cos(2u) = È, , 
È 2 4 

ove si è posto y=0. ‘Infine introducendo la funzione 


ba gin S 


otteniamo un'equazione per la funzione U 


(81) 


(82) 


8 
1 4eV 
+ ( 23"? n? cos(2u))U = 0 


Questa è un'equazione di Mathieu nell'incognita U con autovalore 


8È 


iagt 


1 
np? 


Per il significato della funzione U è necessario che essa risulti periodica di 


periodo 7, mentre per il fatto che È è una correzione al primo ordine in e del- 


l'autovalore, richiediamo che È sia nell'ordine di eV. 


La condizione di periodicità della funzione U può essere soddisfat- 


ta se-e solo se oppure 


(83) 


x = a m= 0,1,2,..., 


2m 


nel qual caso U è la funzione di Mathieu pari 


(84) 


oppure se 


(85) 


U(u) ni cen nin) ms O.l:2si 


\ = bom Mm. 1,2 


nel qual caso U è la funzione di Mathieu dispari 


Per il significato degli autovalori a b,, che compaiono in (83) ed (85) ri- 


2m° “2m 


mandiamo a [7]. 
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Esaminiamo i] caso in cui valgono (83) ed (84) (il caso della funzio 


ne dispari è del tutto analogo); la condizione che È v eV fornisce una condi- 


1 

zione di quantizzazione sulla m della forma 

Vel 5 VeV , 3 
(86) ST = <m<4 an + 8 
cui corrisponde l'autovalore [6] 

r.(m.e,h) 
A eV eV + 
(87) x = 8 E: 7 + 16 ci 7 RE MR i 1 + 0(1) per h>o 
p PURE 9 r(1/4) 


hp 


Evidentemente (82) ed (87) forniscono (70) mentre la funzione di che compare 
in (67) è data da 


essendo m dato da (86). 
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